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Abstract 

Using the low energy effective Hamiltonian with the generalized factorization, we 
calculate the new physics contributions to i? — > iT'^7r~ , Ktt and Kij' in the Topcolor- 
assisted-Technicolor(TC2) model, and compare the results with the available data. 
By using F^'^iO) = 0.20 ± 0.04 preferred by the CLEO data of B ^ tt+tt" decay, 
we find that the new physics enhancements to B ^ Kr]' decays are significant 
in size, ~ 50% with respect to the standard model predictions, insensitive to the 
variations of input parameters and hence provide a simple and plausible new physics 
interpretation for the observed unexpectedly large B Kr]' decay rates. 
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In two-body charmless hadronic B meson decays, new physics beyond the standard 
model (SM) may manifest itself through large enhancements to those penguin-dominated 
decay modes: decays which are expected to be rare in the SM are found to have large 
branching ratios. These potential deviations may be induced by the virtual effects of new 
physics through penguin and/or box diagrams |jl|, ^. 

In the framework of the SM, the two-body charmless hadronic decays B ^ hih2 [ where 
hi and /i2 are the light pseudo-scalar (P) and/or vector(V) mesons ] have been studied 
systematically by many authors |Q |^, |^, H, . On the experimental side, fourteen B^^ 
PP, PV decay channels have been observed by CLEO, BaBar and BELLE Collaboration 



ig, |n],[T2|,|T3|,|I3: 

B TT^n^, Kn, Kt]', p^n^, p°7r±, un^, K*r], K^n^, (f)K^ . (1) 

By comparing the theoretical predictions with experimental measurements one finds 
the following main points: 

• The effective Hamiltonian with generalized factorization approach generally works 
well to interpret the observed pattern of branching ratios. The penguin effects are 



clearly observed |jT5|. 



There may exist a problem to accommodate the data of tttt and Kir simultaneously p6| 
The r]'K puzzle: the B Krj' decay rates are much larger than what one ordinarily 



expected in the SM |T^, [T 



Since 1997, the unexpectedly large branching ratio of i? ^ Krj' has stimulated much 



interests in literature fl^, |T8|, |T9|, ^ . In order to accommodate the data, one may need 
an additional contribution unique to the rj' meson in the framework of the SM, or new 
physics enhancements from new physics models beyond the SM to explain the B — > Krj' 



puzzle |jT2|. In a previous paper |2^, we considered the second possibility and calculated 



the new physics effects on the branching ratios and CP-violating asymmetries of fifty 
seven B PP,PV decay modes in the Topcolor-assited Technicolor (TC2) modelp3[], 
and found that the new physics enhancement to the penguin-dominated decay modes 
can be significant. In another paper we calculated the new physics contributions to 



branching ratios of seventy six B hih2 decay modes in the general two-Higgs-doublet 
models (models I, II and III). 

In this letter, we concentrate on the new physics contributions to seven observed B 
decay modes: B -k^-k^^Kti and B —>■ Krj' in the TC2 model. Particular attention 
is devoted to the details of B ^ Ktt and B Krj' decays when a smaller Fq^'^(O) = 
0.20 ± 0.04 instead of the ordinary F(f ^(0) = 0.33 are used. 

The effective Hamiltonian for the two-body charmless decays B — > hih2 are now known 
at next-to- leading order (NLO) [B^, Rl 0]. The standard theoretical frame to calculate the 
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inclusive three-body decays b sqq^ is based on the effective Hamiltonian[| 



eff 



(AE = 1) = ^ I {v^tV:,Q] + V^tKlQ'j) - Vt,V:, 



10 



J=3 



(2) 

where the operator basis contains the current- current operators Qi,2, the QCD penguin 
operators Q3-6, the electroweak penguin operators Qg-io and the chromo-magnetic dipole 
operator Qg, the exphcit expressions can be found easily for example in Ref.p]. Following 
Ref. 0, we also neglect the effects of the electromagnetic penguin operator Qy^, and do 
not consider the effect of the weak annihilation and exchange diagrams. The coefficients 
Ci in Eq.(|^) are the well-known Wilson coefficients. Within the SM and at scale Mw, the 
Wilson coefficients Ci{Mw) , ■ ■ ■ ^ Cio{Mw) at NLO level and Cg{Mw) at LO level have 
been given for example in pB| . 

Following the same procedure as in the SM, it is straightforward to calculate the new 
7-, Z^- and gluonic penguin diagrams induced by the exchanges of unit-charged scalars, 
the technipion vrf ,7r^ and top-pion tt^ appeared in the TC2 model 0. 

After including the new physics (NP) contributions induced by new penguin diagrams, 
the Wilson coefficients Ci{Mw) i = 1, ■ ■ ■ , 10 at NLO level and Cg at leading-order (LO) 
can be written as 
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(3) 
(4) 

(5) 
(6) 
(7) 
(8) 

(9) 
(10) 



^For b — > dqq decays, one simply make the replacement s ^ d. 

^For more details of TC2 model|^ and the corresponding constraint on its parameter space from the 
data, one can see Refs.[p3|, p6| and our previous paper For the sake of simplicity, we here do not 
present the details about the calculations of new penguin diagrams in the TC2 model, but the reader can 



find them in Ref.|22, M 



3 



sin^ Ow 



105o(xO-4Co(xO+4C^ 



ATP 



Cg{Mw] 



iNP 



(11) 

(12) 
(13) 



where Xt = m'^/M^, the functions -80(2^), Co{x), Do{x), Eq{x) and E'q are the famihar 
Inami-Lim functions which describe the contributions from the VT-penguin and Box di- 
agrams in the SM and can be found, for example, in Refs. p5|. 



The functions Cq^, 

D^P, E^P and E'^^ describe the new physics contributions to Wilson coefficients in the 



TC2 model as given in Ref.p2 
r^NP _ 1 
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TjNP 
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where yt = m* /m\ with m* 
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The first term in Eg .(1141) arises from the top-pion penguins, while the second and third 
term correspond to the color-singlet and color-octect technipion penguin respectively. For 
all four functions, the top-pion penguins always dominate absolutely|2^ . 

In numerical calculations, we use the following parameters of the TC2 model as input 
parameter. Since the new physics contributions from technipions vrf and vr^ are much 
smaller than those from top-pion vr^ within the reasonable parameter space, we here fix 

200GeV for the sake of simplicity: 



m^rj = lOOGeV and m^rg 



mOGeV, 



200GeK F^ = 50GeV, F^ = 120^6^ e = 0.05, 



(24) 
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where and are the decay constants for technipions and top-pions, respectively. 
For m^-, we consider the range of mj^ = 200 ± 100 GeV to check the mass dependence 
of branching ratios of two-body charmless hadronic B meson decays studied. All other 
relevant input parameters, such as the quark masses and form factors, etc., are given in 
the Appendix. 

Since the heavy charged pseudo-scalars appeared in TC2 model have been integrated 
out at the scale Mw, the QCD running of the Wilson coefficients Ci{M]y) down to the 
scale n = 0(7711,) after including the NP contributions will be the same as in the SM. In 
the NDR scheme, by using the input parameters as given in Appendix and Eq.(p^) and 
setting /i = 2.5 GeV, we find that: 

Ci = 1.1245, C2 = -0.2662, C3 = 0.0195, C4 = -0.0441, 

C5 = 0.0111, Ce = -0.0535, Cj = 0.0026, Cg = 0.0018, 

Cg = -0.0175, Cw = 0.0049, Cf = 0.3735 (25) 

where Cf = Cg + C5. 

In this letter, the generalized factorization ansatz^ as being used in Ref. |P, ^ will be 
employed. For the studied seven B meson decay modes, we use the decay amphtudes as 
given in Ref. [§] without further discussion about details. We focus on estimating the new 
physics effects on these seven measured decay modes. 

In the NDR scheme and for SU{3)c, the effective Wilson coefficients can be written 

as m 



eff 




f y + Iv log 



C, + ^ a: (a + C, + Cg) + '^ B[C, ,(26) 



where = (0,0,-1,3,-1,3,0,0,0,0)^, = (0,0,0,0,0,0,1,0,1,0)^, the matrices fy 
and 7y contain the process-independent contributions from the vertex diagrams, and can 
be found, for example, in Refs.f^, |29[. The function Ct, Cp, and Cg describe the contribu- 
tions arising from the penguin diagrams of the current- current Qi,2, the QCD operators 
Qs-Qe, and the tree-level diagram of the magnetic dipole operator Qsg, respectively. 
The explicit expressions of the functions Ct, Cp, and Cg can be found for example in 
Refs.[|, |2^. We here also follow the procedure of Ref. [0] to include the contribution of 
magnetic gluon penguin. 

In the generalized factorization ansatz, the effective Wilson coefficients Cf^ will appear 
in the decay amplitudes in the combinations, 

a2.-i = + a,,^C^f + ^, (z = l,...,5) (27) 

where the effective number of colors N^^ is treated as a free parameter varying in the range 
of 2 < A^^''^ < 00, in order to model the non-factorizable contribution to the hadronic 



^For recent discussions about the generalized factorization approach, one can see Ref.||2^ and reference 
therein. 
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matrix elements. Although A^^''^ can in principle vary from channel to channel, but in the 
energetic two-body hadronic B meson decays, it is expected to be process insensitive as 
supported by the data 0. As argued in Ref.|T7|, Nf{LL) induced by the {V - A){V - A) 
operators can be rather different from N^^{LR) generated by {V — A){V + A) operators. 
In this paper, however, we will simply assume that Nf^{LL) = N^^{LR) = N^^ and 
consider the variation of N^^ in the range of 2 < N^^ < oo since we here focus on the 
calculation of new physics effects. 

In the B rest frame, the branching ratios B{B —>■ PP) can be written as 

where Ttot{B-) = 3.982 x IQ-^^ GeV and 1^4(5°) = 4.252 x IQ-^^GeV obtained by using 
t{B~) = 1.653ps and T(i?°) = 1.548ps |^0|, Pb is the four-momentum of the B meson, 
Mb = 5.279 GeV is the mass of Bu or Bd meson, and 



= 2M^V " + My)2][M| - {Mx - My)2] (29) 

is the magnitude of momentum of particle X and Y in the B rest frame. 

For the seven studied B meson decay modes, currently available measurements from 



CLEO, BaBar and Belle Collaborations O, 12, O, nl are as follows: 



(30) 



+ / (4.3;i|± 0.5) X 10-6 [CLEO], 

B(R K^n^) - / (11-6^2:?1l3) X 10-6 [CLEO], 

- I (l8.8ti^±2.3)xl0-6 [Belle], ^^^^ 

f (17.2111 ±1.2) X 10-6 [CLEO], 
B{B^K+TT-) = I (12.5l|[! 1|:?± 2.3) X 10-6 [BaBar], (32) 

[ (17.4l|^± 3.4) X 10-6 [BELLE], 

B{B^K%+) = (18.2^^0 ± 1-6) X 10-6 [CLEO], (33) 
r^(n^ kKo. S (14.611? 11^) X 10-6 [CLEO], 

^ ^ = [ (2ll?:il||) X 10-6 [BELLE] ^^^^ 

B(B^K^v') = i (80-9° ± 7) X 10-6 [CLEO], 
^ " \ (62 ±18 ±8) X 10-6 [BaBar], 



(2ll?:il||)xl0-6 [BELLE], 

(35) 

B{B^K^ri') = (891}^ ± 9) X 10-6 [CLEO]. (36) 



The measurements of CLEO, BaBar and BELLE Collaboration are consistent with each 
other within errors. 

In Table 1, we present the theoretical predictions of the branching ratios for the seven 
B — s> PP decay modes in the framework of the SM and TC2 model by using the form 
factors from Baner, Stech and Wirbel (BSW) model[^ and Lattice QCD/QCD sum rule 
(LQQSR) model |3^ form factors, as listed in the first and second entries respectively. 
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Table 1: Branching ratios (in units of 10 ^) of seven studied B decay modes in the SM 
and TC2 model by using the BSW and LQQSR form factors, with fc^ = ml/2, A = 0.81, 
A = 0.2205, p = 0.12, r] = 0.34, Nf = 2, 3, oo and = 200 GeV, and by employing 
generalized factorization approach. 
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Theoretical predictions are made by using the central values of input parameters as given 
in Eq.(p4D and the Appendix, and assuming rriji- = 200GeV and A^^''^ = 2, 3, oo in the 
generalized factorization approach. The branching ratios collected in the tables are the 
averages of the branching ratios of B and anti--B decays. The ratio 63 describes the new 
physics correction on the decay ratio and is defined as 



B{B XYf^^ - B{B XY)^^' 

B{B xvy 



6B{B XY) = . ^^.^SM (37) 



By comparing the theoretical predictions with data, the following points can be un- 
derstood: 

• For B^ TT+TT^ decay, the SM prediction is clearly larger than the CLEO mea- 
surement, but agree with BaBar measurement, while the BaBar measurement has a 
larger error than CLEO. The new physics contribution to this tree-dominated decay 
mode is negligibly small. 

• For four B Kn decays, the SM predictions are agree with experimental mea- 
surements. In TC2 models, the theoretical predictions are generally larger than the 
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data but still agree with the data with 2a errors||22| since both the theoretical and 
experimental errors are still large now. 



For —>■ K^Tj' and B^ — > K^rj' decay, the SM predictions are clearly much 
smaller than the data (especially the CLEO measurement). But the new physics 
enhancement can make the theoretical predictions in the TC2 model become agree 
with CLEO/BaBar data within one standard deviation. 



The unexpectedly large B — > Kr]' decay rates were firstly observed in 1997 [Q, 
and confirmed recently with the full CLEO II and II. V samples [|1^. The earlier SM 
predictions in the range of (1 — 2) x 10^^ are too small compared with experiment. In 



the framework of the SM, the B Kr]' decays can be enhanced through [|T6| (i) the 
small running mass at the scale mi, (ii) the sizable SU{3) breaking in the decay 
constant /o and fs, (iii) larger form factor F^^ (0) due to the smaller r] — r]' mixing 
angle —15.4° rather than ^ —20°, (iv) contribution from the rj' charm content, and (v) 
constructive interference in tree amplitudes. However, as pointed out in Ref . ||T8|, 0, the 



above mentioned enhancement is partially washed out by the anomaly effects in the matrix 
element of pseudoscalar densities, an effect overlooked before. As a consequence, the net 
enhancement is not very large: B{B^ K^t]') = (40 — 50) x 10~^ as given in Ref. Ell . 



In the TC2 model, on the other hand, the new gluonic and electroweak penguins 
contribute through constructive interference with their SM counterparts and consequently 
provide the large enhancements, ~ 50% with respect to the SM predictions, as shown in 
Table 0. By using F^'^(0) — 0.33 and other input parameters as given in the Appendix, 
one finds numerically that 

f.(r.± ^ r.±n j (20 - 52) X 10-6 in SM, 

^ ^ ^ ) = \ (30 - 71) X 10-6 in TC2, ^^^^ 

n(RO^K%') - / (19-52) X 10-6 in SM, 

^ ^ ^ ) - \ (29 - 73) X 10-6 in TC2, ^^^^ 

where the effects induced by the uncertainties of major input parameters have been taken 
into account. The SM prediction is still smaller than the CLEO result but agree with 
the BaBar measurement^. In Fig.|l|, we plot the mass- dependence of B{B^ — > K^rj') and 
B{B^ K^i) in the SM and TC2 model. The short-dashed line in Fig.|] shows the SM 
predictions with N'^^ = 3. The dot-dashed and solid curve refers to the branching ratios 
in the TC2 model for A^^''^ = 3 and oo, respectively. The upper dots band corresponds 
to the data with 2a errors: B{B^ K^rj') = (75 ± 20) x 10-6 (average of CLEO and 
BaBar result) and B{B^ K^rj') = {89ttl) x 10-^ (CLEO only). It is evident that the 
theoretical predictions in the TC2 model are agree well with CLEO/BaBar data within 
one standard deviation. 



'^However, a rather small rus is not consistent with recent lattice calculations. 

^One should note that the error of BaBar result is much larger than that of CLEO result. 
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Since —>■ tt+tt^ decay is a tree-dominated decay mode, the new physics correction 
induced through loop diagrams should be very small, as shown in Table |1]. The CLEO 
measurement of this mode puts a very stringent constraint on the form factor F^'^{0): 
F^'^{0) = 0.20 ± 0.04 as given in Ref . [0^ . In the SM, this smaller form factor will lead to 
two difficulties: 

1. First, the predicted B Kir branching ratios will be too small when compared with 
the data since their decay rates depend on the form factors Fq^'^(O) and F^^{0). We 
know that the form factor F^^{0) cannot deviate too much from F^'^(O), otherwise 
the SU (3)-symmetry relation F^'" = Fq'^ will be badly broken. 

2. Second, the predicted B —>■ Krj' branching ratios will be also too small in the SM 
since the branching ratio B{B Krj') depends on both the form factor Fq^'^(O) and 
Fq^ . A small F^^ifi) leads to a small Fq^ and in turn small branching ratios of 
B — > Krj' decays. If we use the relation §], 

^FSi^^.^^), . <r - (40) 

to define Fq"^' with = -9,1° and Os = -22.2° g, the SM prediction for the 
branching ratio B{B —>■ Krj') will be decreased by about 26%. In Table |^, we show 
the branching ratios of seven studied decay modes obtained by using F,f '^(0) = 0.20 
instead of -Fq^'^(O) = 0.33 while keep all other input parameters remain the same as 
being used in Table 0. 

In TC2 model, however, the decrease induced by using smaller F^'^ifd) will be com- 
pensated by large new physics enhancement and therefore restore the agreement between 
the theoretical predictions and the data, as illustrated in Fig.^ for the decay B Krj' . 
By using F^^{<d) = 0.20 ± 0.04, one finds that 

8(B-^KW) = /(12-«)xl0-» in SM, 

^ " \ (20 - 62) X 10-6 in TC2, ^ ' 

r^(RO ^ r^o^,. / (11 - 44) X 10-6 in SM, 

'3iB ^Krj) = I (20 - 61) X 10-6 in TC2, ^^^^ 

where the effects of major uncertainties have been taken into account. 

In Fig.0, we plot the mass dependence of B{B~^ — > K~^rj') and B{B^ — > K^rj') in 
the SM and TC2 model by using F(f^(0) = 0.20 instead of 0.33 (while all other input 
parameters are the same as in Fig.|lD. The short-dashed line in Fig.^ shows the SM 
predictions with A^^''^ = 3. The dot-dashed and solid curve refers to the branching ratios 
in the TC2 model for A^^''^ = 3 and oo, respectively. The upper dots band corresponds 
to the CLEO/BaBar data with 2cr errors. It is easy to see that (a) the gap between the 
SM predictions of B Krj' decay rates and the data is enlarged by using Fq'^{0) = 0.20 
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Table 2: Branching ratios (in units of 10 ^) of seven studied B decay modes in the SM and 
TC2 model by using the BSW form factors with Fo^^(O) = 0.20 instead of Fo^^(O) = 0.33, 
assuming P = ml/2, A = 0.81, A = 0.2205, p = 0.12, r] = 0.34, A^'^^*^ = 2, 3, cx) and 
= 200 GeV, and by employing generalized factorization approach. 
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instead of 0.33, and (b) the new physics enhancement therefore becomes essential for the 
theoretical predictions to be consistent with CLEO/BaBar result within 2cr errors. 

From Table |^, it is easy to see that the new physics enhancement to first three B — > Ktt 
decays and B —>■ Krj' decays are still large in size and play an important role to boost the 
corresponding branching ratios to be consistent with experimental measurements. 

For the B K^tt^ decay mode, the SM predictions are always smaller than the data 
although the error of the data is still very large, as can be seen from Tables 1-2. For 
the case of using F(f''(0) = 0.33, the new physics contribution in TC2 model provide a 
(10 — 20)% enhancement. For the case of using Fq^'^(O) = 0.20, however, the new physics 
contribution in TC2 model result in a (2 — 12)% decrease. We currently are not sure that 
whether there is a discrepancy between the theory and the data for this decay mode. This 
is an open problem now, further refinement of the data will clear this point soon. 

In short, we here studied the new physics contributions to the seven observed B —>■ PP 
decay modes by employing the effective Hamiltonian with generalized factorization. In 
this letter, particular attention is devoted to the details of B ^ Krj' decays, and to the 
discussions about currently known mechanisms to enhance this decay mode. We made the 
numerical calculation by using both F^^{0) = 0.33 as given in the ordinary BSW model, 
as well as F(f^(0) = 0.20 ±0.04 preferred by the CLEO data of B ^ tt+tt-. We presented 
the numerical results in Tables 1-2 and Figs. 1-2. We also discussed the difficulties induced 
by using the smaller F^'^(O) and shown that one can accommodate the data of tt~^tt~, 
K^n, and K^7t~^ simultaneously after taking into account the new physics contributions. 
But we are still not sure if there is a discrepancy between the theory and the data for 
B —>■ K^n^ decay mode. 

By using whether Fq^'^(O) = 0.20 ±0.04 or 0.33, we always found that the new physics 
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enhancements to -B ^ Krj' decays are significant in size, and hence the theoretical pre- 
dictions of B{B Kt]') in the TC2 model are agree with CLEO/BaBar data within 
2(7 errors. This seems to be a simple and plausible new physics interpretation for the 
observed rj'K puzzle. 



Appendix: Input parameters 

In this appendix we present the relevant input parameters. 

• Input parameters of electroweak and strong coupling constant, gauge boson masses, 
light meson masses, ■ ■ ■, are as follows (all masses in unit of GeV )||, 



= 1/128, as{Mz) = 0.118, sin^ = 0.23, Gp = 1.16639 x 10-^(G'eV")-^ 
Mz = 91.188, Mw = 80.42, m^o = m^± = 5.279, m^± = 0.140, 
m^o = 0.135, = 0.547, m^/ = 0.958, mK± = 0.494, m^o = 0.498. (43) 

For the elements of CKM matrix, we use Wolfenstein parametrization, fix the pa- 
rameters A, X and p to their central values: A = 0.81, A = 0.2205, p = 0.12, but 
varying rj in the range of ?7 = 0.34 ± 0.08. 

We first treat the internal quark masses in the loops as constituent masses: m?, = 
4.88GeV, rric = 1.5GeV, rris = 0.5GeV, rriu = ma = 0.2GeV. Second, we use the 
current quark masses for (i = u, d, s, c, b) which appear through the equation of 
motion when working out the hadronic matrix elements. For p = 2.5GeV, one found 
0: rub = 4.88GeV,mc = 1.5GeV,m, = 0.122GeV, ma = 7.6MeV, m„ = 4.2MeV. 
For the mass of heavy top quark we use rrit = rn^{mt) = 168GeV. 

For the decay constants of light mesons, the following values will be used in the 
numerical calculations (in the units of MeV): 

U = 133, fK = 158, /; = = 78, = /J = 68, 

= -0.9, f^;, = -0.23, = -113, f^, = 141. (44) 

where /^(/^ and /^(z, have been defined in the two- angle-mixing formalism with 6q = 



-9.1° and 9s = -22.2= 
The relevant form factors are ||^ 

Fo^-(O) = 0.33, F^^{0) = 0.38, Fo^''(0) = 0.145, Fo^'''(0) = 0.135, (45) 



in the BSW model 311, and 



Fo^-(O) = 0.36, F(f^(0) = 0.41, Fo^''(0) = 0.16, Fq^^^O) = 0.145, (46) 
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in the LQQSR approach. And the momentum dependence of form factor -Fo(^^) 
was defined in Ref.pl| as -Fo(A;^) = Fo(0)/(l — A;^/m^(0+)). The pole masses being 
used to evaluate the fc^ dependence of form factors are m(0+) = 5.73 GeV for ub 
and db currents, and m(O^) = 5.89 GeV for sb currents. 
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Figure 1: Plots of branching ratios of decays 5+ — > K'^r]' (la) and — > K^rj' 
(lb) versus mass in the SM and TC2 model with Fo^''(0) = 0.33 . The 
short-dashed hne shows the SM predictions with Nf^ — 3. The dot-dashed 
and sohd curve refers to the branching ratios in the TC2 model for A^^*'^ = 3 
and oo, respectively. Theoretical uncertainties are not shown here. The dots 
band corresponds to the CLEO/BaBar data with 2a errors: B{B^ K^rj') — 
(75 ± 20) X 10-6 and B{B^ K\) = {89ttl) x lO'^. 
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Figure 2: Same as Fig.|but for F(f^(0) = 0.20 instead of 0.33. 
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